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Low-noise temperature measurements at frequencies in the milli-Hertz range are required in the 
LISA (Laser Interferometer Space Antenna) and LISA PathFinder (LPF) missions. The required 
temperature stabiUty for LISA is around 10 /^K Hz"^''^ at frequencies down to 0.1 mHz. In this paper 
we focus on the identification and reduction of a source of excess noise detected when measuring time- 
varying temperature signals. This is shown to be due to non-idealities in the ADC transfer curve, 
and degrades the measurement by about one order of magnitude in the measurement bandwidth 
when the measured temperature exhibits drifts of ~ /^Ks~^. In a suitable measuring system for the 
LISA mission, this noise needs to be reduced. Two different methods based on the same technique 
have been implemented, both consisting in the addition of dither signals out of band to mitigate the 
ADC non-ideality errors. Excess noise of this nature has been satisfactorily reduced by using these 
methods when measuring temperature ramps up to 10 /iKs"^. 

PACS numbers: 07.20.Dt, 07.87.+V 
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I. INTRODUCTION 

LISA (Laser Interferometer Space Antenna) is a joint ESA-NASA space mission conceived to detect Gravitational 
Waves (GWs) [1]. LISA consists in a constellation of three spacecraft in the vertexes of an equilateral triangle 5 million 
kilometres to the side. The constellation orbits the Sun following the ecliptic, some 20 degrees (45 million kilometres) 
behind the Earth. Each spacecraft houses two proof masses in nominal free fall, and laser links are established between 
spacecraft — see Figure [TJ The role of the links is to enable interferometric measurements of relative distance and 
acceleration variations between pairs of proof masses in distant spacecraft. This is how LISA will detect GWs, since 
they show up as tidal forces (or geodesic deviations in the language of General Relativity Theory) in the region where 
the detector is. 

Expected GW signals are extremely weak [2| , hence severe requirements must be imposed on the proof mass residual 
acceleration noise to ensure LISA meets its GW detection objectives. These are defined in terms of spectral density Q: 
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in the frequency band from 0.1 mHz to 0.1 Hz, which is where optimum GW detection performance can be obtained 
with interferometer arm- lengths of ^--^5 million kilometres [3]. 

The conditions in Equation ([T]) are not only very demanding, they cannot be directly put to test in an earth based 
environment. The reason is of course the impossibility to maintain a physical system in accurate free fall near the 
Earth's surface during periods of several hours, as required for measurements in sub-milli-Hcrtz frequency bands. 
This fact motivated ESA, the European Space Agency, to fly a technology precursor mission to secure proper working 
of the key technologies needed by LISA. Such technology mission goes by the name LISA PathFinder (LPF), and its 
launch date is scheduled for 2011. LPF is a squeezed version of one LISA arm: its length is downscaled to ~30cm, 
and the two proof masses are housed in a single spacecraft. Noise requirements in LPF are relaxed by about an 
order of magnitude relative to those of LISA. Also, the Measuring Bandwidth (MBW) in the LTP (LISA Technology 
Package, the main instrument on board LPF) is reduced to the range between 1 mHz and 30 mHz. Because of its small 
dimensions, LPF cannot work as_a milli-Hertz GW detector; rather, LPF is a noise monitor intended to understand 
and properly model its sources ,5;]. 

One of such sources of noise, which can limit the performance of LISA, is temperature fluctuations. Thermal 
stability is necessary to ensure the stability of the optical elements properties in the Optical Metrology System 
(Optical Bench, mirrors, beam-splitters, etc.), as well as of the proof masses' environment, where such effects as 
radiometer and/or radiation pressure fluctuations will cause random accelerations if the temperature fluctuates. This 
happens in LISA and LPF alike, and the maximum tolerable levels of temperature fluctuations can be derived from 



2 




estimates of the magnitude of their contribution to the total instrument noise, and (conventionaUy) assuming that 
they should contribute not more than ^10% to the total. This generates a requirement in temperature stability for 
LISA of [H 

4^lisaM<10/^KHz-i/2 (1) 

in LISA'S MBW, i.e., from 0.1 mHz to 0.1 Hz. Consequently, measurements capable to discern such small temperature 
fluctuations are needed. In LPF, the requirements being relaxed as described above, the temperature stability is also 
less severe Q: 

(w) < 100 /iK Hz"i/2 (2) 

in the frequency band from 1 mHz to 30 mHz. In order to make meaningful temperature measurements, the measuring 
system must be quieter than this figure, and again 10% of the above is required for it, i.e., 

4(systcmH<10MKHz-l/2 (3) 

The LTP temperature monitoring system has already been designed and tested 01 — see also Section|lI]below — and 
it works. However, a small excess noise has been detected at frequencies around the milli-Hertz when the measured 
temperature drifts with time. This effect does not pose a serious problem in the LTP, as it only shows up at very 
low frequencies (< 1 mHz) and in the presence of somewhat high temperature drifts, \dT/dt\^ /iKs^^. In LISA the 
MBW stretches down to 0.1 mHz and even weaker temperature drifts are likely to significantly deteriorate the 
performance of the measuring system. As we shall see, its origin has been traced to non-linearities of the Analog-to- 
Digital Converters (ADC). 

This paper focuses on the identification and reduction of this excess noise. It is organised as follows: Section [H] 
briefiy describes the temperature measurement system which will actually fiy in LPF. Section IIIII delves into the 
details of a potential problem with the non-idealities of the ADCs. In Section IIVI we present experimental evidence 
that the observed excess noise can indeed be attributed to the ADCs non-ideal response, and also two different ways 
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to reduce that noise are analysed. In Section |V] we give details of the experimental setups, tests carried through to 
assess the quality of the proposed methods, and results. Finally, Section IVTl highlights the most relevant results and 
the prospects of their applicability for LISA. 



II. TEMPERATURE MEASUREMENT SYSTEM DESCRIPTION 

In this section the temperature measurement system is briefly reviewed. Details can be found in ^] . The functional 
block diagram of the system is given in Fig. [21 
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FIG. 2: Temperature measurement system chain. 



The temperature signal, x{t), coming from a Wheatstone bridge is modulated by a square wave signal of angular 
frequency ujc {uJc/2tt — 6.25IIz) to avoid the 1// noise of the Instrumentation Amplifier (lA). The modulated signal 
is amplified by the lA and low-pass filtered to avoid aliasing. The signal is then quantised by a 16-bit ADC. Once 
the signal is digitised, a digital demodulation takes place: N^v samples are averaged during either polarity of the 
modulating square wave, and the averaged value of one polarity is subtracted from that of the other one — see 
schematics in Fig. [2j This process results in a noise equivalent temperature for dc signals of 

at frequencies down to 0.1 mHz 0, [I]- 

The signal coming from the analog circuit of the front-end electronics (FEE) is quantised by an ADC. Only one 
type of ADC suitable for this purpose has been certified for use in space applications; it is a 16-bit capacitor based 
Successive Approximation Register (SAR). Limitations in the performance of the system at low frequencies (milli- 
Hertz range) when slowly drifting input signals (~ s^^) are quantised are related to the non-idcalitics of such 
ADCs, as will be shown below. 



III. NON-IDEAL QUANTISATION NOISE 



In this section we analyse the ADC non-linearity errors, and how they perturb the temperature measurements. 

Quantisation is inherently a non-linear process. By construction, the values of the analog function are rounded up 
to the closest ADC step, so that even an ideal ADC generates output errors associated to the differences between the 
real and quantised values of the signal. In an ADC the step size, or Least Significant Bit (LSB), is defined by 

^-^Iv^^FS (5) 

where N is the number of bits of the ADC, and Vps is the maximum voltage (full-scale) the ADC can quantise. 
When a large number of bits is considered the error introduced by the quantisation process is usually treated as an 
independent random variable with uniform probability density function (pdf ) and white spectral density Q . 

In a real ADC, however, the quantisation steps are not uniform due to mismatches in the internal topology of the 
ADC; more specifically, tolerances in the capacitors of the SAR array — see Figure [3] Such non-uniformity 
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is specified with two parameters: Differential Non-Linearity (DNL) and Integral Non-Linearity (INL) errors. The DNL 
is defined as the deviation between two adjacent transition points of the quantisation ladder and an LSB, or (2^. [23| 

DNL{i)^K,+i-K,-A, z = 0,1, 2,. ..,2^-1 (6) 

where Ki is the output value of the i-th quantisation code. 

The additional noise related to the DNL errors can be reduced by suitable dithering [l^ [ij, [l^ [l^ . The inherent 
ADC noise plus FEE noise can be considered in this case as a dither source, which appears to be enough to make the 
DNL negligible in practice — see Section fill CI below. 

On the other hand, the INL error is defined as the discrete integral of the DNL, i.e., 

i 

INL{i) = Y,DNL{j), 0,1, 2,. ..,2^-1 (7) 

i=o 

and can be understood as the difference between the real and the ideal ADC transfer curves — see Figure [3] The noise 
introduced by this error is usually less noticeable, although much more difficult to reduce. In this paper we focus 
on how to do this, since the INL is the one limiting the performance of the measurement when slowly drifting input 
signals are present. This is a common situation in LTP temperature measurements. 
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FIG. 3: INL error effect on a 8-bit ADC: the INL is the difference between dots (ideal ADC) and crosses (real ADC). 



Experimental results — see Section fill Dl — obtained with the temperature measurement system reflect the problem 
depicted in Figure [31 which shows up as excess noise in the spectrum. Such noise could not be attributed to ac- 
tual temperature fluctuations since sensors were placed in a thermal environment where fluctuations were efficiently 
screened out Jj], hence such behaviour may not be due to variable temperatures. Analog noise or interferences coming 
from the signal processing chain (Wheatstone bridge, amplification, low-pas filtering) were also discarded as possible 
sources of the extra noise. Finally, the excess noise disappears when a dc temperature value (i.e., not drifting with 
time) is measured [l^ — see Figure [51 All in all, the origin of the excess noise has been identified as due to INL errors 
of the ADC. In the following we present more quantitative arguments confirming this hypothesis. 



A. ADC bit error description 



We first analyse the effect of a faulty bit in a SAR ADC on the performance of the measurement. For this, a simple 
ADC model is considered. The analog-to-digital conversion is done (in SAR ADCs) by comparing the sampled signal 
with an analog voltage generated by a Digital-to- Analog Converter (DAC) and a SAR [10]. The topology of the DAC 
is based on a switching capacitor bank composed by 16 capacitors, scaled from 2^^C to C, where 2^^C defines the 
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Most Significant Bit (MSB) and C defines the LSB. The DAC output voltage is ^ 



VDAC = ''^FS — —FTTn TTTT ^FS 



(2 



DC 



(8) 



where bk is the binary digit {k = Q stands for the LSB and fc = — 1 for the MSB); it is set to or 1 depending 
on the sampled voltage, and N is the number of bits of the ADC. An error 5Ck = Cfc, real— 2*^(7 in the fc-th capacitor 
results in an error in the DAC output when the corresponding bit is set to 1. The voltage error efc induced by a faulty 
k-th bit is therefore 



bk SCk - . , SCk . 



(9) 



where A is the ideal LSB. 

The erroneous bit produces a superimposed periodic pattern in the quantisation error of the ADC with 1 LSB of 
amplitude when the input voltage varies between and Vpsi and exhibits a periodicity with the ADC input voltage 
of 2*^"'"^ A. An example of this is shown in Figure S] where it can be seen how the faulty bit introduces a long period 
component superimposed on the typical sawtooth error function of an ideal ADC. 
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FIG. 4: Three simulated 4- bit ADCs. Top: quantisation error functions for the ideal ADC (black trace), a real ADC with a 
faulty k = l bit (ei = 0.3A, red), and a real ADC with a faulty fe = 3 bit (ei = 0.6A, blue). Centre and bottom: error differences 
between the two real ADCs and the ideal one, respectively. Units in ordinates are LSBs. 



B. Dither signal effect in ideal quantisers 

When a dither voltage, d, with a certain probability density function (pdf), p{d), is added to the input signal of 
the ADC, V, the average quantisation error observed at the ADC output, {q{v)), is given by 



ili^)) = / q{v + z)p{z)dz 

J —00 



(10) 
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where q(v) is the quantisation error of an ideal ADC. If we define the Fourier transform in voltage domain by [2a | 

/•OO 

Q(0= / q{v)e-'^^dv (11) 



and apply it to Ea. ((TO| . the following ensues: 

(Q(0) = S(0^*(0 (12) 

The quantisation error for an ideal ADC (a sawtooth waveform of amplitude A/2 — see Figure S]) is a periodic 
function of v with period A. It can be expanded in Fourier series: 

, , A ^ (-1)"+! . 2Tmv 
TT '^-^ n A 

n— 1 



and its Fourier transform is 

n — — OO 
(n#0) 

where 5{-- ■) is Dirac's 5-function. If we consider zero-mean Gaussian noise as the dither voltage, its pdf and its 
Fourier transform are, respectively, 

^ vki^'^ '''' 
no = (16) 

where cr^ is the dither variance. Substituting now Eqs. and into Eq. and taking the modulus, we obtain 



1(2(6)1 = A-e-^^^^/^ E -^f^-^) (17) 



71 

From Eq. ([T7]) we note that the Gaussian dither signal, d, low-pass filters the quantisation error — see Figure [51 top. 
For instance, for = A the attenuation of the first term in the series is 2.7x 10~^. In our case, we have a 16-bit ADC 
with Vfs=10V, gwasi-white noise at the input of the ADC of ~ 7/xVHz-i/2 and a noise-equivalent bandwidth 
(NEBW) of 1.2 /cut-off = 600 Hz, where /cut-off is the cut-off frequency of the anti-ahas filter Thus, A and a are 
readily calculated, i.e., 

A = ^^ = 0.15mV (18) 

a = sI/^{lo) X (1.2 /cut-off)'/' = 0.17 mV (19) 

In our specific case, cr ~ A, hence the quantisation noise from the ideal ADC should be suppressed by the inherent 
noise of the analog processing chain which acts as a filter to the ideal quantisation error. This can be very clearly 
seen in Figure [5l top, which shows the errors of an ideal ADC (vertical lines) alongside the equivalent low-pass filter 
profile (dashed trace) achieved by means of Gaussian dither with a ~ A. 



C. Dither signal effect in non-ideal quantisers 

The error of a real ADC is formed by the ideal quantisation error function — see Eq. ([13]) — plus a term related to 
the non- idealities of the ADC — see Figure [H Thus, 

q{v) ^ qi{v) + qk{v) (20) 

where qi(v) and qk{v) are the ideal quantisation error and the quantisation error due to the faulty bits, respectively. 
As seen in the previous section, the ideal quantisation error is filtered out by the analog noise in the measurement 
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chain, which acts as a dither signal of Gaussian pdf. Instead, the quantisation error due to the non-ideaHty of the 
ADC is not reduced by the same analog noise, which causes it to show up as extra noise. In this section we show how 
it can be identified in the temperature measurements, and how it limits the performance of the system. 
The modulus of the Fourier transform of the non-ideal noise for a single defective bit, k, say, is [16| 
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where is given in Eq. ([9]) . The effect of Gaussian dither on qk (v) is readily calculated with Eq. ([12 
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Equations (pij) and are plotted in Figure [5] (centre and bottom) . The plot in the centre corresponds to errors 
associated to fc = 0, i.e., the LSB. Here, the fundamental period is 2A. The plot in the bottom shows the same for 
fc = 3 (the 4-th bit) with a fundamental period of 16A. For each plot, the low-pass filter generated by Gaussian 
dither of (T = A is also plotted (dashed lines). It can be noticed that this dither signal suffices to suppress the noise 
associated to a faulty LSB, fc = (this corresponds to the DNL effect which, is readily mitigated). However, it cannot 
attenuate the low frequency lines due to errors in the higher bits. In fact, when fc = 3 the dithering with cr = A is 
clearly insufficient. 
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FIG. 5: Top: quantisation error for an ideal ADC, QiiX)- Centre: quantisation error for a real ADC with fc = as the faulty 
bit. The error is eo = 0.5A. Bottom: same as previous but for fe = 3. The error is e3 = 0.5A. Note that the x-axis scale of the 
bottom plot is different from the top and middle ones; for instance, the main frequency for Qo is 1/2A, while that for Qa is 
1/2'' A. This means that errors in the higher bits show up at lower frequencies than those due to errors in the LSB. 



From Figure [5] we note how faulty bits introduce spectral components in the measurement if the signal spans a 
large enough fraction of the range of the ADC transfer curve. Conversely, if the input signal is a constant dc then no 
extra noise will be seen, since no bit, whether faulty or not, will change its state. Therefore, the INL effect becomes 
perceptible in our system only when the input signal runs through a sufficiently wide fraction of the ADC range. 



D. INL effect on general signals 

As we have just seen, when the input signal, w, is not a dc constant (e.g., a ramp), the errors related to the INL of 
the ADC tend to introduce low frequency noise components which degrade the performance of the system. In fact, 
this is seen to happen in our device for voltage variations in the order of few milli-volts. Moreover, if the input signal 
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changes rate, i.e., i' 7^ then the frequency components introduced by the INL errors spread across the frequency 
band. 

In order to deal with non-constant signals and to estimate INL induced noise, a simplification is expedient. First, 
let us assume the input signal is a straight line for a certain time interval, i.e., 

( a + bt if < i < t,nax 

vit) = { 23) 

[0 if i < or t > imax 

where a is the initial value of the signal v(t\ and 6 is the slope of the input signal (Vs""'^). If we note qhiyi) the usual 
time-frequency Fourier transform of the non-ideal quantisation noise, it immediately follows from Eq. (|23p that 

iwa/fc { u:\ 



In this case, Eq. ([2T|1 becomes 



n— — 00 I ^ ^ ^ 

and the main frequency component for an ADC with an error in the A;-th bit is located at 

Before we proceed further, a technical comment is in order. In several equations above, Dirac ^-functions appear. 
They are the result of infinite length integration intervals in Fourier transform calculations, which are of course 
idealisations. In any practical situation, such intervals are limited to the experimentally available data ranges, so that 
the (5-functions are actually sinc-functions: they have identical centre points but spread around those centres depending 
on the integration interval lengths. In order to make sense of e.g. Eq. ()25|) we must ask which is the minimum required 
integration time to obtain a meaningful spectrum or, in other words, which is the minimum duration, tmax of the 
ramp signal in Eq. (j23p . This is easily inferred from the frequency of the lowest spectral line uj\^ in Eq. ([26]): for a 
conventional ten cycle integration time we get 

27r 2*^+^ A 

imax > 10 = 10 — — (27) 

These considerations apply to the analysis of the ADC response to signals which drift slowly with time, where slow 
drift means such signals can be conveniently approximated by a series of concatenated ramps with suitable slopes, 
and lengths complying with Eq. ([27|) . In these circumstances, we can generalise Eq. ([26]) as follows: 



uJi,k = ^ \v{t)\ , tl^i.fc = ^ m (28) 
From Eqs. and (pS)) some conclusions can be drawn: 

• High-slope signals at the ADC input translate into high-frequency components at the output of the ADC. 

• Errors in the higher bits show up as noise at low frequency. 

• Consequently, high-slope input signals combined with errors in high bits may show as noise peaks in the mea- 
surement bandwidth. 

• The fundamental frequency associated to a faulty bit varies with the variations of v{t), hence the error in a 
faulty bit spreads across the frequency band when v ^ 0. 

The above can be validated by looking at real data from our LTP temperature measurement system. We took a long 
time series of temperature data, about 4x10^ seconds, and subdivided it into shorter stretches of ~16000 seconds. 
With this, we calculated the Short-Time Fourier Transform (STFT) [isj . also known as spectrogram, and plotted the 
results as shown in Figure [6] — see caption for details. The data are amenable to the analysis described above, and 
the results indicate the presence of the foreseen frequency peaks associated to faulty bits. They are clearly visible in 
FigurelH bottom graph, as the darker areas, which show how the fundamental frequency of the INL errors for different 
faulty bits is varying with time, and accurately tracing the \v\ curve shown on the top left plot. A straightforward 
fit to the first equation Ij28p permits the identification of the corresponding faulty bit, thus confirming that the extra 
noise in the temperature measurements is indeed due to the INL errors of the ADC. 

Once the problem has been identified, we focus on possible solutions to fix it. The following section describes two 
different methods tested successfully to avoid the consequences of INL errors in real ADCs. 
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FIG. 6: Top left: Absolute value of the input signal slope, [^(t)!, of the measurement. It varies between and 12.5 /iVs^^. The 
sensitivity at the input of the ADC is ~1.35VK~^, hence the values of the voltage slope almost directly provide the actual 
temperature drifts in [Ks~^]. Top right; power spectral density for: (i) the signal shown in the top left plot, i.e., a non-rfc 
signal (solid black trace) and, (ii) a dc signal (dashed red trace), i.e., with \v{t)\ < 0.5 /iVs~^. The excess noise when measuring 
non-constant signals is very noticeable, and degrades the performance of the measurement by about one order of magnitude. 
Bottom: STFT (or spectrogram) of the measurement. The energy of the signal is concentrated in specific frequencies which 
change with time precisely following |i'(t)|, as predicted by Ea. (|28|l . The latter is represented by dashed superimposed traces, 
and are labelled by the order of the corresponding faulty bit (from k = 2 to k = 7). Experimental results and theoretical 
estimates are in excellent agreement. 

IV. MITIGATION OF INL ERRORS 

Two different techniques have been used to deal with the non-linearities of the ADC. The first is based on the 
injection of Gaussian noise out of the MBW, taking advantage of the oversampling involved in the measurement. The 
second is based on the addition of a triangular wave of high frequency to the signal quantised by the ADC [23| . The 
following sections describe both techniques as well as their practical implementation. 

A. Gaussian noise signal injection 

In Section IIIIBI we saw that Gaussian noise can be used as a dither signal in order to mitigate the non-idealities 
of the ADC transfer curve. We have also seen that the amount of noise generated by the measurement system itself 
(cr~A) is not enough to suppress the periodic components of erroneous bits for k>2 — see Fig. [6] Thus, a natural 
solution is to add more Gaussian noise to the ADC input. Obviously, this noise should be added out of the MBW in 
order not to disturb the frequency range of interest, which is in the milli-Hertz range. The ADC sampling frequency 
is 38.4kHz, thus leaving a large frequency slot to accommodate the required additional noise. Care must however be 
taken to also place it away from the fundamental frequency of the modulating square wave (and its harmonics) so as 
to avoid bringing the Gaussian noise back into the MBW in the demodulation process. The required amount of noise 
to be injected, characterised by ct, basically depends on the input signal slope, v(t), and on the frequencies of interest, 
and is limited by the digital processing performed after quantisation — see below. Figure [7] shows the relationship 
between the input signal slope and the main frequency component for each of the faulty bits. It is useful to identify 
the faulty bits potentially affecting the measurement as a function of the input signal drifts. 
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FIG. 7: Relationship between the input signal slope, the faulty bit and the fundamental frequency affected by the error in the 
k-ih bit. When high input signal slopes are present, more faulty bits will affect the measurement quality in the MBW. Clearly, 
LSBs show up at higher frequencies than MSBs. The shaded areas span the MBWs of LISA (green) and LPF (yellow), as 
indicated. Note that LPF's MBW is a sub-region of EISA's. 



Once the problematic bits are identified, the needed amount of Gaussian dither must be calculated. For this, we 
define the a needed to filter out the noise components associated to the A;-th faulty bit. If we (conventionally) adopt 
a damping factor of 10 for the first harmonic then a is easily derived from the condition (see Eq. ([H])) 
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Gaussian noise dither: practical implementation 



In this section we briefly describe the hardware implementation of the Gaussian dithering scheme. The circuit 
diagram is given in Fig. [S] the first stage amplifies the noise of an operational amplifier (OP-07 [IJl); the second 
stage is a high-pass filter of 4-th order (two Sallen-Key filters in cascade) with a cut-off frequency of 100 Hz; the third 
stage is an adder that sums the signal of interest (the amplified output of the Wheatstone bridge) to the dither signal. 
In the fourth stage, the sum of signal and dither are low-pass filtered with a 4-th order low-pass filter (again, two 
Sallen-Key filters in cascade) with cut-off frequency 3 kHz. The output of this chain is fed to a 16-bit SAR ADC. The 
noise shape of the dither signal is given in Fig. [9] (solid trace) where the characteristic frequencies of the high- and 
low-pass filters are clearly visible. 

As mentioned in Section |TT1 the amount of Gaussian dither must be kept under control to avoid excessive noise 
folding back into the MBW during the digital demodulation stage. The limit of the usable Gaussian noise amplitude, 
cr, can be easily estimated assuming a flat spectrum, S'y, dither, from Whi/Ztt (the corner frequency of the high-pass 
filter) to lom/^tt (the corner frequency of the low-pass filter). The demodulation by a square wave of frequency lOc 
entails a certain gain at dc of the odd harmonics of uJc of the dither noise which, under the just mentioned assumption, 
results in the following p^ : 



<S'v, extra (^^ = 0) ~ S'y^ dither ' 



^ (2^ 

n=[w,-„/2aJcJ 



1) = 



(31) 
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FIG. 8: Gaussian noise dither signal generator. Its power spectral density is shown in Fig. (9] v{t) is the amplified signal from 
the Wheatstone bridge — see Figure [2] d{t) is the dither signal before low-pass filtered. 



where [ J is the floor function, e.g., [x\ is the largest integer smaher than x. 

In order to cut down the noise leaking into the MBW, i.e., at dc in practice, we again impose a requirement that 
it be less than 10% of the floor noise in the absence of ADC errors, S'y, fee, i-e., 



Sv, cxtra(w — 0) < 



10 



(32) 



In our case, (^1-^/271 ~ 100 Hz, lom/^tt ~ 3 kHz, and lOc/2t: = 6.25 Hz, so that the S-sum in Eq. (|3ip is approximately 
0.015, hence = 0) ~ 0.015 5'y, dither- The value of S'v.fee is ^ 1.5 x lO-^V^Hz-i in the MBW. The maximum 

value for 5*^, dither compliant with the criterion Eq. ([32|) is therefore 10~^" Hz~^, which corresponds to cr ~ 0.5 mV. 
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FIG. 9: Solid red trace: Noise shape of the dither signal generated by the circuit shown in Figure (8) The central plateau 
corresponds to cr = 3 mV — see Section IV Bl The dashed black trace is the floor noise of the measurement system itself, prior 
to demodulation Q. 



In Section lYl we present the tests and results obtained with this technique. 



B. Triangular wave dither 



Another implemented and tested technique consists in using a deterministic signal instead of random Gaussian noise 
as the dither signal. Different signals can be used for this purpose, but for simplicity we have used a triangular wave 
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form in the generation of the signal — see below. In this section we describe the theoretical basis and the hardware 
implementation details. 

The dither signal is now a triangular wave, which can be represented by the Fourier series 



1 



(2n+ 1)2 



cos {2n + l)ujtit 



(33) 



where Do and Wtr are the amplitude and the frequency of the triangular wave, respectively — see Fig. [TT] Its pdf 
is M 



p{d) 



d < and d > Do 



(34) 



whose Fourier transform is 



sm{Do/2)^ 



As shown in Section [lll|, the averaged quantisation error is 

(Q(?)> = Q(C)^*(0 

Eq. ([36|) is graphically evaluated in Figure [10] for different faulty bits and triangular wave amplitudes. 
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FIG. 10: Top: quantisation error for a real ADC with a faulty bit in fc=0, eo = 0.5A, and the low-pass filter when using a 
triangular wave as the dither signal with Do = 2A. Centre: same as above for an error in bit fc=3, ea = 0.5A. Bottom: same 
as above but with Do = 8 A instead of 2 A. Note that a;-axis scales are different for Qo and Q3. 



Like we did with the Gaussian noise, we can give an expression to determine the required triangular wave amplitude, 
Do^ to attenuate the effect of the fc-th bit: 



sin (L'o/2)f 



(^o/2)e 



1 27r 
< — for f = -r-mr 
-10 ^ 2'=+! A 



(37) 



There is no closed form solution to Eq. (|37p for the amplitude of the triangular wave, which requires numerical 
evaluation in each specific case. The amplitude Do needed to attenuate the errors coming from the fc=6 bit is ~20 mV. 
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Triangular wave dither: practical implementation 

Figure [11] shows the triangular wave added to the system. It is important to note that the addition of the triangular 
signal will not perturb the temperature measurement within the MBW. The digital demodulation involved in the 
measurement — see Sec. |Tl] — is done by averaging 6144 samples: 3072 during one polarity and 3072 during the 
opposite polarity. Afterwards they are subtracted and divided by 2. Thus, if we inject exactly the same signal in both 
polarities the net contribution of the added dither signal to the output is zero — see Eq. ([55)) . 



Vpos(t) 




FIG. 11: Triangular wave dither signal and signal coming from the measurement chain (dashed trace). 
The output signal after the digital processing is 

[t^posft) + d{t)] - Kcgit) + d{t)] Vpos - Vncg ,„„^ 

= ^ = ^ (38) 

where d{t) is the triangular wave — see Fig. [11] — and an overbar (~) means average over 3072 samples. 

Eq. (|38p indicates that there is no limit on the amplitude of the triangular dither, provided triangular waves in both 
polarities are identical. Therefore that non-linearities of MSBs can be reduced without degrading the measurement. 
This happens for an analog dither signal, but ours is actually generated with a DAC, which imposes some limits on 
the validity of the previous statement — see below. 

The circuit which implements the triangular signal is shown in Figure 1121 and consists in an 8-bit up-and-down 
counter followed by a 12-bit DAC where only the 8 LSBs are used. The circuit is configured such that the quantisation 
step of the DAC closely matches that of the 16-bit ADC. Mismatches here result in actual performance diverging 
from the predictions of the theoretical analysis described above. 

The output signal from the DAC is added to the signal of interest and both are low-pass filtered (a 4-th order 
Sallen-Key filter with a cut-off frequency of 3kHz), then quantised by the 16-bit SAR ADC. The triangular wave is 
low-pass filtered to eliminate high-frequency components related to the digital quantisation of the DAC. Nevertheless, 
the previous analysis is still valid since the fundamental frequency of the triangular wave is 50 Hz and the low-pass 
filter cut-off frequency is 3 kHz. Thus, the dither signal will go almost unaltered through the filter (except for some 
distortion in its high-frequency components). 




CLK 



FIG. 12: Triangular wave dither signal generator. v{t) is the amplified output signal of the Wheatstone bridge. d{t) is the 
triangular wave. 
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V. TEST SET-UP AND RESULTS 



Both methods described in Section IIVI have been put to test. For comparison, a 24-bit Dcha-Sigma ADC 
(LTC2440 which, in principle, should exhibit less non-linearity problems, has also been tested. The two dither 
techniques under study have been tested using exactly the same electronics, the same 16-bit SAR ADC AD977 [HI [13, 
and under the same input signal conditions. In this section first we give a brief description of the test set-up and, 
then, the obtained results are presented. 



A. Test set-up 



The test set-up is composed by different parts. Temperature sensors are placed inside a thermal insulator designed 
to screen out ambient temperature fluctuations to the required level in the MBW, i.e., S'^'^^(a;) < 10/iKHz~^/^ for 
a;/27r > 1 niHz A temperature control is included to implement different temperature profiles inside the insulator; 
basically, a set of temperature ramps with different slopes is generated in order to assess whether or not the methods to 
overcome the INL errors of the ADC work — see Figure fT4l The temperature control consists in a heater commanded 
by a programmable power supply which is in turn controlled by the computer calculated value of the difference between 
the desired temperature and the actual measurement — see Figure [T3l 




thermal insulator temperature measurement system 



FIG. 13: Setup scheme. FEE labels the analog signal processing of the temperature measurement system and PS stands for 
the programmable power supply. 



The nominal and actual profiles are shown in Figure [Ml The experiment has been repeated for the two dithering 
techniques, and for the 24-bit Delta-Sigma ADC as well. Prior to that, the 16-bit SAR ADC with no dither signal 
was tested in order to provide a reference measurement with the INL effects in it. 




FIG. 14; Top: nominal and experimental (achieved with the feedback temperature control shown in Fig. I13p temperature 
profiles. Bottom: time derivative of the signals shown in the top plot. Tested slopes are: 0.5, 1, 2, 4, 8 and 16/iKs~^. 
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B. Experimental results 



The power spectral density in the different tested configurations (no dither, Gaussian noise, triangular wave and 
24-bit Delta-Sigma ADC) are given in Fig.[Tn]for three different slopes, more specifically, for 1, 4 and 8/iKs~^. The 
measurements performed with the 16-bit SAR ADC are clearly affected by the INL errors when no dither signal is 
used. The noise in the MBW increases by more than one order of magnitude when slopes are around 1 /zKs~^ and 
above. Figure [T5l shows how the INL errors appear in the MBW, spanning wider frequency regions as slopes increase, 
thus confirming the predicted behaviour. For instance, when the slope is 4/iKs~^ the INL effect is only noticeable at 
frequencies below 6mHz, while for a slope of 8/LtKs^^ the noise appears at frequencies as high as lOmHz. 

When the dither signal is Gaussian noise, the INL effects can be satisfactorily reduced down to 1 mHz for drifts 
under 8/iKs^^. Looking up Figure [7l we see that faulty bits up to fc = 5 will create additional noise in the band. 
Then, using Eq. ([30]) . we find that the necessary dither requires cr~ 22 A, or (T~3mV. This is however a factor of 6 
larger than the maximum estimated after Eq. (|32p . which means some extra noise will be added to the system floor 
noise, if we insist on applying a cr = 3 mV dither, as discussed in Section [IV Al We did take this option, with the result 
that the floor noise becomes a factor ~1.5 larger than nominal [in good agreement with predictions calculated using 
Eq. ([3T|) ]. with the advantage that good damping of the ADC's INL errors obtains. 
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FIG. 15: Linear Power Spectral Density for different configurations: 16-bit ADC without dither, 16-bit ADC with Gaussian 
dither, 16-bit ADC with triangular wave dither, and 24-bit Sigma-Delta ADC without dither. Input signals are ramps with 
slopes varying from 1 to 8 /^Ks"^. In absence of dither, the noise of the ADC increases (in amplitude and bandwidth) with the 
slope of the input signal. The solutions adopted work satisfactorily for slopes up to ~8^Ks~^ level. The 24-bit Sigma-Delta 
does not exhibit non-linearity problems for the tested slopes. See text in Section [V Bl for details. 



The use of a triangular wave as the dither signal appears as the most robust option to deal with the non-idealities 
of the ADC: on the one hand, the floor noise is left untouched and, on the other hand, immunity to faulty bits 
and high signal slopes can be tuned essentially at will. The 24-bit Delta-Sigma ADC exhibits a behaviour similar 
to that observed in the measurements performed by the 16-bit ADC with added dither signals. Non-linearity errors 
in this ADC are not noticeable when measuring signals drifting up to 8/LtKs~^. The generated triangular wave had 
the following properties: -Do=155mV and a period of 20 ms, which is an integer sub-multiple of the duration of each 
polarity (80 ms) — see Figure [TT] This amplitude is almost 8 times higher than the one estimated towards the end 
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of Section flV Bl (20 mV) to attenuate errors up the k = 6 bit. This means that immunity to faulty bits under higher 
temperature drifts is accomphshed. Figure \W\ provides a clear display of the superiority of triangular wave over 
Gaussian noise dither. 
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FIG. 16: Gain of the equivalent filters of the Gaussian noise dither and the triangle wave. The gain shown for each bit 
corresponds to the fundamental frequency. Gaussian noise dither mitigates the error in the bits fc < 5 whereas the triangular 
wave attenuates the error in the bits fc < 9. 



Summing up, Figure [TSl confirms the analysis in this paper, and shows that when slopes of ^lO^Ks""'^ are present, 
the effects of the non-idealities can be made negligible by use of proper dither. For higher slopes, e.g., 16/iKs~^, 
some increase in the power spectrum near one milli-Hertz is detected. We believe this is probably due to the relative 
simplicity of the model used. Further research is ongoing to clarify these matters. Nevertheless, such drifts look 
unlikely both in LISA and LPF. 



VI. DISCUSSION 



Temperature diagnostics measurements in LISA and LISA PathFinder must meet very demanding requirements: 
very low noise, a few /iKHz"^/^, and very low frequency band, below 1 Hz and down to fractions of a milli-Hz. After 
a suitable measurement system was in place, extra noise was seen to appear at the lower end of the MBW when slowly 
drifting temperatures were measured. For temperature slopes higher than ~0.5 /iKHz~^/^ this extra noise challenges 
the performance of the LTP measuring system, the problem growing more severe at lower frequencies, down into the 
submilli-Hertz LISA band. 

The source of the noise has been investigated and identified as due to the INL effects of the ADC, which is a 16-bit 
S AR ADC [1^ . We have laid down the theoretical basis of the problem, which has been validated by laboratory 
experiment. 

Once the problem is well understood, it needs to be solved. Two different options to mitigate this effect have been 
proposed, analysed and tested. Both are based on the addition of a dither signal to the signal of interest prior to 
the ADC quantisation. In one of them, we use Gaussian noise, which certainly reduces the INL error effects of the 
ADC, although in our case its performance is limited by the the digital processing which sets an upper limit on the 
noise amplitude, a, which can be injected in the measurement. Alternatively, dithering with a triangular wave does 
not add noise in the MBW, and generally results in a more efficient attenuation of INL errors. A 24-bit Delta-Sigma 
ADC has been used for comparison, with the result that a 16-bit ADC plus averaging (to increase the resolution) and 
proper dither reaches the same performance as a 24-bit Delta-Sigma ADC. 

Summing up, a robust method to suppress the INL errors of the ADC in the temperature measurement system 
of the LPF mission has been presented in this paper. Such a system is not foreseen to fly in LPF, but the results 
obtained herein are important in view of the LISA mission which will need similar (or even higher) precisions at lower 
frequencies. The method we have described is also useful for other subsystems of LISA where ADC non-linearities 
appear as a limitation in the performance of the measurement system. 
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